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FOREWORD 


This  final  report  is  submitted  in  fulfillment  of  CDRL  CLIN  0001,  Data  Item  A001,  Title: 
Scientific  and  Technical  Reports  of  a  Small  business  Technology  Transfer  (STTR)  Phase  II 
contract  No.  F33615-96-C-3217  entitled,  “Development  of  the  Aerodynamic/Aeroservoelastic 
Modules  in  ASTROS,”  covering  the  performance  period  from  24  September  1996  to  24 
September  1998. 

This  work  is  the  second  phase  of  a  continuing  two-phase  STTR  contract  supported  by 
AFRL/Wright-Patterson.  The  first  phase  STTR  contract  No.  F33615-95-C-3219  entitled, 
“Enhancement  of  the  Aeroservoelastic  Capability  in  ASTROS,”  was  completed  in  May  1996 
and  published  as  WL-TR-96-31 19. 


Both  STTR  Phase  I  and  Phase  II  contracts  are  performed  by  the  same  ZONA  Team  in  which 
ZONA  Technology,  Inc.  is  the  prime  contractor,  whereby  the  team  members  include:  the 
University  of  Oklahoma  (OU),  Universal  Analytics,  Inc.  (UAI),  and  Technion  (I.T.T.). 

This  final  report  consists  of  eight  volumes,  these  are: 


ASTROS* 
Volume  I 
Volume  II 
Volume  in 
Volume  IV 


ZAERO  User’s  Manual 
ZAERO  Programmer’s  Manual 
ZAERO  Application  Manual 
ZAERO  Theoretical  Manual 


ASTROServo 
Volume  I 
Volume  II 
Volume  IE 
Volume  IV 


Aeroservoelastic  Discipline  in  ASTROS,  User’s  Manual 
Aeroservoelastic  Discipline  in  ASTROS,  Programmer’s  Manual 
Aeroservoelastic  Discipline  in  ASTROS,  Application  Manual 
Aeroservoelastic  Discipline  in  ASTROS,  Theoretical  Manual 


This  document  (Volume  II)  is  the  Programmer’s  Manual  of  the  Aeroservoelastic  (ASE) 
interaction  module  developed  to  facilitate  ASE  analysis  and  the  application  of  ASE  stability 
and  response  constraints  within  ASTROS. 


At  AFRL/Wright-Patterson,  Captain  Gerald  Andersen  was  the  contract  monitor  and  Dr.  V.  B. 
Venkayya  was  the  initiator  of  the  whole  STTR  effort.  The  technical  advice  and  assistance 
received  from  Mr.  Doug  Niell  of  the  MacNeal  Schwendler  Corporation,  Dr.  V.  B.  Venkayya 
and  others  from  AFRL  during  the  course  of  the  present  phase  on  the  development  of 
ASTROS*  are  gratefully  acknowledged. 
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Chapter  1 


Introduction 


The  Programmer's  Manual  describes  the  new  engineering  application  modules  and  data  base 
entities  used  by  the  new  aeroservoelastic  (ASE)  module. 
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Chapter  2 


Engineering  Application  Modules 


2.1  Summary  of  new  application  modules 

The  engineering  application  modules  added  to  support  the  ASE  module  are: 


ASEDRV 

ASEFLUT 

ASERUN 

ASESENS 

BDMAT 

CNTMOD 

GMAT 

GMGKUPD 

GSENS* 

GSTMOD 

MGFL 

MIST 

RMODE 

UPDV 


driver  for  ASE  options 

ASE  analysis  w/  or  w/o  control  effects 

ASE  run  options 

ASE  constraint  sensitivity 

baseline  dynamic  matrices 

aerodynamic  forces  and  mass  matrices  associated  with  control 

baseline  structural  matrices 

update  generalized  matrices 

sensitivity  of  generalized  matrices 

aerodynamic  forces  associated  with  gust  modes 

update  dynamic  matrices 

aerodynamic  approximations 

reduced-size  modal  analysis 

design  variable  update 


*  -  a  MAPOL  procedure 


A  detailed  description  of  the  application  modules  is  given  in  Section  2.2. 


modes 
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2.2  Description  of  application  modules 

The  engineering  application  modules  added  to  support  the  ASE  module  are  listed  in  Table  2.1. 


Module 

Baseline 

Analysis 

Restart 

Optim. 

Cold 

Optim. 

Description 

ASEDRV 

X 

X 

X 

X 

driver  for  ASE  options 

ASEFLUT 

X 

X 

X 

X 

flutter  w/  or  w/o  control  effects 

ASERUN 

X 

X 

X 

X 

ASE  run  options 

ASESENS 

X 

X 

flutter  constraint  sensitivity 

BDMAT 

X 

X 

baseline  dynamic  matrices 

CNTMOD 

X 

X 

control  modes  matrices 

GMAT 

X 

(X) 

X 

baseline  structural  matrices 

GMGKUPD 

X 

X 

X 

update  generalized  matrices 

GSENS* 

X 

X 

sensitivity  of  gen.  matrices 

GSTMOD 

X 

X 

gust  modes  matrices 

MGFL 

X 

X 

X 

update  dynamic  matrices 

MIST 

X 

(X) 

X 

aerodynamic  approximations 

RMODE 

X 

X 

X 

reduced-size  modal  analysis 

UPDV 

X 

X 

design  variable  update 

(X)  -  only  in  final  analysis 
*  -  a  MAPOL  procedure 


Table  2.1:  Engineering  application  modules  in  the  ASE  module 


3 


Engineering  Application  Module:  ASEDRV 
Entry  Point:  ASEDRV 


Purpose:  MAPOL  director  for  aeroservoelastic  analysis 
MAPOL  calling  sequence: 

CALL  ASEDRV  (BC,  SUB,  LASE,); 

BC  Boundary  condition  number  (Integer,  Input) 

SUB  Flutter  subcase  number  to  be  processed  in  this  path  (Integers,  Input) 

LASE  If  TRUE,  perform  ASE  discipline.  Otherwise  perform  FLUTTER  discipline 
(Logical,  Output) 

Application  calling  sequence: 

None 

Method: 

The  FLCOND  ID  for  the  SUB  flutter  case  is  extracted  from  the  CASE  relation  and  used  to 
extract  the  METHOD  string  from  the  associated  FLUTTER  template.  If  METHOD=ASE, 
LASE  is  set  to  TRUE.  Otherwise,  LASE  is  set  to  FALSE  and  the  other  output  parameters  are 
irrelevant. 

Design  Requirements: 

1.  This  module  is  the  driver  for  all  modules  of  the  ASE  discipline. 

Error  Conditions: 

1.  None. 


4 


Engineering  Application  Module:  ASERUN 
Entry  Point:  ASERUN 
Purpose:  Driver  for  ASE  options. 

MAPOL  calling  sequence: 

CALL  ASERUN  (ASEDIS,  ASEDBMAK,  ASEOPTDB,  NUP); 


ASEDIS 

Logical  flag  indicating  that  ASE  discipline  is  to  be  employed 
(Logical,  Output) 

ASEDBMAK 

Logical  flag  indicating  that  a  modal  data  base  is  to  be  created  in  an 
analysis  run  (Logical,  Output) 

ASEOPTDB 

Logical  flag  indicating  that  a  modal  data  base  is  to  be  created  in  the 
first  iteration  of  the  current  optimization  run  (Logical,  Output) 

NUP 

Number  of  design  iterations  between  modal  data  base  updates 
(Integer,  Output) 

Application  calling  sequence: 

None 

Method: 

The  module  reads  from  relation  ASEUP  the  values  of  the  parameters  ASEDB  and  NUP.  If 
the  number  of  entries  in  relation  ASEUP  is  zero,  the  logical  flag  ASEDIS  is  set  to  FALSE 
to  tell  the  MAPOL  sequence  that  the  ASE  discipline  is  not  employed.  If  ASEDB=1,  logical 
flag  ASEDBMAK  is  set  to  TRUE.  If  ASEDB=0,  logical  flag  ASEOPTDB  is  set  to  TRUE. 

Design  Requirements: 

1.  None. 

Error  Conditions: 

1 .  Relation  ASEUP  contains  more  than  one  entry. 
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Engineering  Application  Module:  ASEFLUT 


Entry  Point:  ASEFLUT 
Purpose: 

To  perform  ASE  analysis  and  to  calculate  ASE  constraints  in  optimization  run  with  applied 
ASE  constraints. 


MAPQL  calling  sequence: 


CALL  ASEFLUT  (  NITER,  BC,  SUB,  ESIZE(BC),  GSIZEB,  NRSET,  MHHFL(BC,SUB)], 
[KRHHFL(BC,SUB)],  [BHHFL(BC,SUB)],  [EMS(BC,SUB)], 
[DTMS(BC,SUB)],  [AMS(BC,SUB)],  [RMS(BC,SUB)], 

[MHC(BC,SUB)],  [PfflG(BC)];  [NLCR(BC,SUB)],  [PHFLO(BC,SUB)], 
[PLROG(BC,SUB)],  ASELAMBD,  ASEMARG,  CONST ); 


NITER 

Design  iteration  number  (Integer,  Input) 

BC 

Boundary  conditions  identification  number  (Integer,  Input) 

SUB 

Flutter  subcase  number  (Integer,  Input) 

ESIZE(BC) 

The  number  of  extra  point  degrees  of  freedom  (Integer,  Input) 

GSIZEB 

The  size  of  the  structural  set  (Integer,  Input) 

NRSET 

The  number  of  supported  degrees  of  freedom  (Integer,  Input) 

[MHHFL(BC,SUB)] 

Modal  mass  matrix  (Input) 

[KRHHFL(BC,SUB)] 

Real  modal  flutter  stiffness  matrix  (Input) 

[BHHFL(BC,SUB)] 

Modal  flutter  damping  matrix  (Input) 

[EMS(BC,SUB)] 

Matrix  E  of  aerodynamic  approximation  (Input) 

[DTMS  (B  C,  SUB)] 

Matrix  DT  of  aerodynamic  approximation  (Input) 

[AMS  (BC, SUB)] 

Matrix  which  consists  of  matrices  A\  of  aerodynamic 
approximation  (Input) 

[RMS(BC,SUB)] 

Vector  of  lag  terms  of  aerodynamic  approximation  (Input) 
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[MHC(BC,SUB)] 

Mass  coupling  matrix  for  structural  and  control  modes  (Input) 

[PfflG(BC)] 

Eigenvalues  in  the  global  set  (Input) 

[NLCR(BC,SUB)] 

Vector  of  column  ranges  for  Roger  approximation  (Input) 

[PHFLO(BC,SUB)] 

Vector  containing  the  numbers  of  modes  to  be  used  in  ASE 
analysis  (Input) 

[PLROG(BC,SUB)] 

Truncation  vector  for  Roger  aerodynamic  approximation 
(Output) 

ASELAMBD 

Relation  containing  the  results  of  flutter  analysis  (Output) 

ASEMARG 

Relation  containing  the  results  of  control  margin  analysis 
(Output) 

CONST 

Relation  of  constraints  values  (Output) 

Application  calling  sequence: 

None 

Method: 

The  module  retrieves  the  flutter  discipline  entries  from  relation  CASE  for  the  current  boundary 
conditions.  Then  it  reads  data  from  relations  FLUTTER  and  FLFACT  for  the  current  subcase 
identification  number.  If  the  set  of  omitted  modes  is  not  empty  truncation  of  the  modal  and 
aerodynamic  matrices  is  performed.  Then  two  loops  are  initialized:  outer  loop  for  the  density 
ratios  and  inner  -  for  velocities.  Inside  the  inner  loop  for  each  pair  density-velocity  the  hidden 
module  AEMODL  is  called  to  create  matrix  [^4ae]  of  the  state  space  aeroelastic  model  (Eq.  4.5) 
of  the  Theoretical  Manual).  If  closed-loop  flutter  analysis  is  required  subroutines  PMODEL 
and  MKVEHI  are  called  to  build  plant  and  vehicle  state-space  models,  and  finally  matrix  [Av] 
of  the  closed-loop  system  is  created. 

Eigenvalues  of  the  matrix  [. Aae ]  (open-loop  flutter  analysis)  or  [Av]  (closed-loop  flutter 
analysis)  are  calculated,  sorted  and  finally  output  data  are  written  into  relation  ASELAMBD.  If 
margin  analysis  is  required  subroutines  MARGNS  and  MIMOMR  are  called  to  calculate  SISO 
gain  and  phase  margins,  and  MIMO  margins  which  written  into  relation  ASEMARG.  If  ASE 
gust  analysis  is  required  subroutine  CGUSTR  are  called  to  calculate  gust  response.  If  the 
flutter  analysis  is  performed  in  optimization  run  with  applied  flutter  constraints,  data  from 
relation  DCONFLT  are  retrieved,  the  constraint  values  are  calculated  and  written  into  relation 
CONST.  The  constraint  identification  numbers  are  written  into  relation  ASELAMBD.  If  the 
flutter  analysis  is  performed  in  closed-loop  optimization  run  with  applied  margins  constraints, 
data  from  relation  DCONUGM,  DCONLGM,  DCONUPM,  DCONLPM,  are  retrieved,  the 
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constraint  values  are  calculated  and  written  into  relation  CONST.  The  constraint 
identification  numbers  are  written  into  relation  ASEMARG.  Finally,  the  results  of  the  flutter 
analysis  are  written  to  the  output  file. 

Design  Requirements: 

1.  None. 

Error  Conditions: 

1.  None. 
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Engineering  Application  Module:  ASESENS 

Entry  Point:  ASESENS 

Purpose: 

To  compute  the  sensitivities  of  active  ASE  constraints  in  the  current  boundary  condition. 


MAPOL  calling  sequence: 


CALL  ASESENS  (NITER,  BC,  SUB,  ESIZE(BC),  GSIZEB,  NDV,  [MHHFL(B C, SUB)] , 

[KRHHFL(BC,SUB)],  [BHHFL(BC,SUB)],  [DGKV(BC)],  [DGMV(BC)], 
[DGMVC(BC,SUB)],  [EMS(BC,SUB)],  PTMS(BC,SUB)], 
[AMS(BC,SUB)],  [RMS(BC,SUB)],  [MHC(BC,SUB)], 

[PfflG(BC)],  [PHFLO (B C, SUB)] ,  [PLROG(BC,SUB)], 

ASELAMBD,  ASEMARG,  CONST,  [AMAT] ); 


NITER 

BC 

SUB 

ESIZE(BC) 

GSIZEB 

NDV 

[MHHFL(BC,SUB)] 
[KRHHFL(BC,SUB)] 
[BHHFL(B  C,  SUB)] 
[DGKV] 

[DGMV] 

[EMS  (B  C,  SUB)] 
[DTMS(BC,SUB)] 


Design  iteration  number  (Integer,  Input) 

Boundary  conditions  identification  number  (Integer,  Input) 

Flutter  subcase  number  (Integer,  Input) 

The  number  of  extra  point  degrees  of  freedom  (Integer,  Input) 

The  size  of  the  structural  set  (Integer,  Input) 

Number  of  design  variables  (Integer,  Input) 

Modal  mass  matrix  (Input) 

Real  modal  flutter  stiffness  matrix  (Input) 

Modal  flutter  damping  matrix  (Input) 

A  matrix  of  design  sensitivity  of  generalized  stiffness  matrix 
(Input) 

A  matrix  of  design  sensitivity  of  generalized  mass  matrix  (Input) 
Matrix  [£]  of  aerodynamic  approximation  (Input) 

Matrix  [Z)T]  of  aerodynamic  approximation  (Input) 
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[AMS(BC,SUB)] 

Matrix  which  consists  of  matrices  [A  {\  of  aerodynamic 
approximation  (Input) 

[RMS(BC,SUB)] 

Vector  of  lag  terms  of  aerodynamic  approximation  (Input) 

[MHC(BC,SUB)] 

Mass  coupling  matrix  for  structural  and  control  modes  (Input) 

[PfflG(BC)] 

Eigenvalues  in  the  global  set  (Input) 

[PHFLO(B  C,  SUB)] 

Vector  containing  the  numbers  of  modes  to  be  used  in  flutter 
analysis  (Input) 

[PLROG(BC,SUB)] 

Truncation  vector  for  Roger  aerodynamic  approximation 
(Output) 

ASELAMBD 

Relation  containing  the  results  of  flutter  analysis  (Output) 

ASEMARG 

Relation  containing  the  results  of  control  margin  analysis 
(Output) 

CONST 

Relation  of  constraints  values  (Output) 

[AMAT] 

Matrix  containing  the  sensitivities  of  constraints  to  design 
variables  (Output) 

Application  calling  sequence: 

None 

Method: 

The  module  first  reads  input  matrices  from  the  data  base.  If  the  set  of  omitted  modes  is  not 
empty  truncation  of  the  modal,  aerodynamic  and  sensitivity  matrices  is  performed. 

Then  a  loop  is  initialized  to  retrieve  flutter  constraint  value  and  its  identification  number  from 
relation  CONST.  According  to  the  constraint  identification  number,  the  corresponding  entry  in 
relation  ASELAMBD  is  extracted,  velocity  and  dynamic  pressure  are  determined  and  the 
hidden  module  AEMODL  is  called  to  create  matrix  [Aae]  of  the  state  space  aeroelastic  model. 

If  open-loop  flutter  optimization  is  performed  then  the  hidden  module  DAPDV  is  called  to 
calculate  the  sensitivity  of  the  matrix  [Aae]  to  design  variables,  and  sensitivities  of  the 
eigenvalue  and  the  current  flutter  constraint  are  then  calculated.  If  closed-loop  flutter 
optimization  is  performed  subroutines  AEMODL,  PMODEL,  MKVEHI,  DAVBDV, 
DAVBDG  are  called  to  calculate  sensitivities  of  the  state  space  matrices  of  closed-loop  system 
to  structural  and  gain  design  variables.  Than  sensitivity  of  the  constraint  is  calculated 
according  to  its  type  :  flutter,  control  margin,  or  open-loop  gain.  Finally,  the  sensitivity  of  the 
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constraint  is  added  as  a  column  in  matrix  [AMAT],  the  loop  is  closed  and  the  next  constraint  is 
treated. 


1.  None. 


Error  Conditions: 
1.  None. 
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Engineering  Application  Module:  BDMAT 
Entry  Point:  BDMAT 


Purpose: 

To  calculate  the  baseline  generalized  direct-input  mass,  stiffness,  and  damping  matrices,  and  a 
real  modal  stiffness  matrix  and  for  flutter  analysis. 


MAPOL  calling  sequence: 

CALL  BDMAT  (  NITER,  BC,  SUB,  HSIZE(BC),  ESIZE(BC), 

[GMB(BC)],  [GKB(BC)],  [MHHFL(B C, SUB)] ,  [KHHFL(B C, SUB)] , 
[BHHFL(BC,SUB)],  [KRHHFL(BC,SUB)],  [GM2(BC,SUB)], 
[GK2(BC,SUB)],  [GB2(BC,SUB)] ); 


NITER 

BC 

SUB 

HSIZE(BC) 
ESIZE(BC) 
[GMB(BC)j 
[GKB(BC)] 
[MHHFL(B  C,  SUB)] 
[KHHFL(B  C,  SUB)] 
[BHHFL(B  C,  SUB)] 
[KRHHFL(B  C,  SUB)] 
[GM2(BC,SUB)] 
[GK2(BC,SUB)] 
[GB2(BC,SUB)] 


Design  iteration  number  (Integer,  Input) 

Boundary  conditions  identification  number  (Integer,  Input) 
Flutter  subcase  number  (Integer,  Input) 

The  number  of  eigenvectors  computed  (Integer,  Input) 

The  number  of  extra  point  degrees  of  freedom  (Integer,  Input) 
Generalized  mass  matrix  of  the  baseline  (Input) 

Generalized  stiffness  matrix  of  the  baseline  (Input) 

Modal  mass  matrix  (Input) 

Modal  flutter  (real  or  complex)  stiffness  matrix  (Input) 

Modal  flutter  damping  matrix  (Input) 

Real  modal  flutter  stiffness  matrix  (Output) 

Generalized  direct-input  mass  matrix  (Input) 

Generalized  direct-input  stiffness  matrix  (Input) 

Generalized  direct-input  damping  matrix  (Input) 
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Application  calling  sequence: 


None 

Method: 

The  module  reads  the  modal  flutter  stiffness  matrix  from  the  data  base.  If  the  matrix  is  real,  it 
is  written  to  the  data  base  as  [KRHHFL(BC,SUB)].  If  it  is  complex  its  real  part  is  written  as 
[KRHHFL(BC,SUB)].  Then  module  reads  the  names  of  K2PP,  M2PP,  and  B2PP  for  the 
current  flutter  subcase  from  relation  CASE.  If  the  K2PP  name  is  not  empty,  generalized  direct- 
input  stiffness  matrix  [GK2(BC,SUB)]  is  computed  as  the  difference  between  the  modal  flutter 
stiffness  matrix  [KRHHFL(BC,SUB)]  and  the  generalized  stiffness  [GKB(BC)]  matrix  of  the 
baseline  (Eq.  (2.24)  of  the  Theoretical  Manual).  If  the  M2PP  name  is  not  empty,  generalized 
direct-input  mass  matrix  [GM2(BC,SUB)]  is  computed  as  the  difference  between  the  modal 
flutter  mass  matrix  [MBHFL(BC,SUB)]  and  the  generalized  mass  matrix  [GMB(BC)]  of  the 
baseline  (Eq.  (2.25)  of  the  Theoretical  Manual).  If  the  B2PP  name  is  not  empty  the  module 
reads  the  set  identification  number  of  user-input  damping  parameters  (namely  entities 
VSDAMP  and  TABDMP1)  for  the  current  flutter  subcase  from  relation  CASE.  If  it  is  zero 
direct-input  damping  matrix  [GB2(BC,SUB)]  is  equal  to  the  modal  flutter  damping  matrix 
[BHHFL(BC,SUB)].  Otherwise  user-input  damping  matrix,  which  is  the  sum  of  the  first  and 
the  third  terms  in  the  Eq.  (2.26)  of  the  Theoretical  Manual,  is  created.  Then  direct-input 
damping  matrix  [GB2(BC,SUB)]  is  computed  as  the  difference  between  the  modal  flutter 
damping  matrix  [MHHFL(BC,SUB)j  and  the  user-input  damping  matrix. 

Design  Requirements: 

1.  None. 

Error  Conditions: 


1.  None. 
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Engineering  Application  Module:  CNTMOD 
Entry  Point;  CNTMOD 


Purpose: 

To  calculate  generalized  unsteady  aerodynamic  control  forces  and  control  modes  mass 
matrices. 


MAPOL  calling  sequence: 


CALL  CNTMOD  ( NITER,  BC,  SUB,  [AJC],  [SKJ],  [SCNTLK],  [ACNTLK],  [SCNTLG], 
[ACNTLG],  REUNMK,  [PHKH],  [PffiG(BC)],  [MGG], 
[QHC(BC,SUB)],  [MCC(BC,SUB)],  [MICB(BC,SUB)],  CONTRFL, 
ISYM ); 


NITER 

Design  iteration  number  (Integer,  Input) 

BC 

Boundary  conditions  identification  number  (Integer,  Input) 

SUB 

Flutter  subcase  number  (Integer,  Input) 

[AJC] 

Unsteady  pressure  coefficients  at  the  J-set  due  to  unit  control 
surface  deflections  (Input) 

[SKJ] 

Integration  matrix  relating  the  unsteady  pressure  coefficients  at  the 
J-set  to  the  unsteady  aerodynamic  forces  in  the  K-set  (Input) 

[SCNTLK] 

Displacements  and  slopes  in  the  K-set  due  to  unit  symmetric 
control  surface  deflection  (Input) 

[ACNTLK] 

Displacements  and  slopes  in  the  K-set  due  to  unit  antisymmetric 
control  surface  deflection  (Input) 

[SCNTLG] 

Displacements  and  slopes  in  the  G-set  due  to  unit  symmetric 
control  surface  deflection  (Input) 

[ACNTLG] 

Displacements  and  slopes  in  the  G-set  due  to  unit  antisymmetric 
control  surface  deflection  (Input) 

REUNMK 

Relation  generated  by  the  UZAERO  module  (Input) 

[PHIKH] 

Eigenvectors  in  the  K-set  (Input) 

[PHIG(BC)] 

Eigenvectors  in  the  global  set  (Input) 
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[MGG] 

Mass  matrix  in  the  global  set  (Input) 

[QHC(BC,SUB)] 

Generalized  aerodynamic  control  force  (Output) 

[MCC(BC,SUB)] 

Mass  coupling  matrix  for  control  modes  (Output) 

[MHC(BC,SUB)] 

Mass  coupling  matrix  between  structural  and  control  modes 
(Output) 

CONTRFL 

Logical  flag  indicating  that  control  system  are  employed 
(Logical,  Output) 

ISYM 

Logical  flag  indicating  symmetric  or  antisymmetric  case  is 
employed  (Logical,  Output) 

Application  Calling  Sequence 
None 


Method: 


The  module  reads  flutter  identification  number  for  the  subcase  SUB  from  relation  CASE  and 
then  extracts  control  system  ID  number  from  relation  ASESOL.  If  it  is  zero,  logical  flag 
CONTRFL  is  set  to  FALSE  and  the  module  stops.  Otherwise  logical  flag  CONTRFL  is  set  to 
TRUE  and  matrices  [QHC(BC,SUB)],  [MCC(BC,SUB)],  and  [MICB(BC,SUB)]  are 
calculated. 
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Engineering  Application  Module:  GMAT 

Entry  Point:  GMAT 

Purpose: 

To  generate  the  baseline  diagonal  generalized  stiffness  and  mass  matrices. 
MAPOL  calling  sequence: 


CALL  GMAT  ( NITER,  BC,  LAMBDA,  Mil,  [GMB(BC)],  [GKB(BC)] ); 

NITER  Design  iteration  number  (Integer,  Input) 

BC  Boundary  conditions  identification  number  (Integer,  Input) 

LAMBDA  Relation  of  normal  modes  eigenvalues  output  from  the  REIG  module 
(Input) 


[Mil]  Modal  mass  matrix  created  in  the  REIG  module  (Input) 

[GMB(BC)]  Generalized  mass  matrix  of  the  baseline  (Output) 
[GKB(BC)]  Generalized  stiffness  matrix  of  the  baseline  (Output) 
Application  calling  sequence: 

None 


Method: 

The  module  reads  the  generalized  stiffnesses  from  relation  LAMBDA,  builds  a  diagonal 
generalized  stiffness  matrix,  and  writes  it  into  the  data  base.  Then  it  writes  [Mil]  as  the 
generalized  mass  matrix  of  the  baseline,  [GMB(BC)]. 

Design  Requirements: 

1.  None. 


Error  Conditions: 
1.  None. 
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Engineering  Application  Module:  GMGKUPD 

Entry  Point:  GMGKUPD 

Purpose: 

To  calculate  the  generalized  mass  and  stiffness  matrices  for  the  modified  structure. 


MAPOL  Calling  Sequence: 

CALL  GMGKUPD  (NITER,  NITERUP,  NDV,  DVFLG,  [DGMV(BC)],  [DGKV(BC)], 

[DELDV],  GLBDES,  [GMB(BC)],  [GKB(BC)],  [DELGM],  [DELGK], 
[GM(BC)3,  [GK(BC)j,  [DGMVC(BC,SUB)],  [MICB(BC,SUB)], 
[MIC[BC,SUB)]); 


NITER 

NITERUP 

NDV 

DVFLG 

[DGKV] 

[DGMV] 

[DELDV] 

GLBDES 

[GMB(BC)] 

[GKB(BC)] 

[DELGM] 

[DELGK] 


Design  iteration  number  (Integer,  Input) 

Design  iteration  number  in  which  the  modal  data  base  was  created 
or  updated  (Integer,  Input) 

Number  of  design  variables  (Integer,  Input) 

Logical  flag  indicating  that  design  variables  were  updated  via 
VINIG  bulk  data  entries  (Logical,  Input) 

A  matrix  of  design  sensitivity  of  generalized  stiffness  matrix 
(Input) 

A  matrix  of  design  sensitivity  of  generalized  mass  matrix  (Input) 

A  vector  of  differences  between  updated  design  variables  and 
those  used  for  building  modal  data  base  (Input) 

Relation  of  global  design  variables  (Input) 

Generalized  mass  matrix  of  the  baseline  (Input) 

Generalized  stiffness  matrix  of  the  baseline  (Input) 

A  matrix  of  additions  to  a  generalized  mass  matrix  because  of 
design  variables  updating  (Output) 

A  matrix  of  additions  to  a  generalized  stiffness  matrix  because  of 
design  variables  updating  (Output) 
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[GM(BC)] 

[GK(BC)] 

[DGMVC] 

[MICB  (B  C,  SUB)] 

[MIC(BC,SUB)] 

Application  Calling  Sequence: 

None 
Method: 

The  module  reads  the  baseline  generalized  matrices,  and  their  sensitivities,  from  the  data 
base.  If  NITER  is  zero  or  one  and  the  flag  DVFLG  is  TRUE,  it  also  reads  vector  DELDV. 
If  NITER  >1  the  module  retrieves  the  design- variable  values  from  relation  GLBDES  and 
makes  the  vector  of  differences  DELDV.  Then  the  module  calculates  the  generalized 
matrices  of  the  modified  structure  according  to  the  Eq.  (2.10),  (2.1 1)  of  the  Theoretical 
Manual. 

Design  Requirements: 

1.  None 

Error  Conditions: 


Generalized  mass  matrix  of  the  modified  structure  (Output) 

Generalized  stiffness  matrix  of  the  modified  structure  (Output) 

A  matrix  of  design  sensitivity  of  mass  coupling  between 
structural  and  control  modes  (Input) 

A  matrix  of  mass  coupling  between  structural  and  control  modes 
of  the  baseline  (Input) 

A  matrix  of  mass  coupling  between  structural  and  control  modes 
of  the  modified  structure  (Output) 


1.  None. 
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Engineering  Application  Module:  GSTMOD 
Entry  Point:  GSTMOD 


Purpose: 

To  calculate  generalized  unsteady  aerodynamic  gust  forces. 

MAPOL  calling  sequence: 

CALL  GSTMOD  (NITER,  BC,  SUB,  [QGK],  REUNMK,  [PHEKH],  [QHG(BC,SUB)] ); 
NITER  Design  iteration  number  (Integer,  Input) 

BC  Boundary  conditions  identification  number  (Integer,  Input) 

SUB  Flutter  subcase  number  (Integer,  Input) 

[QGK]  Intermediate  unsteady  gust  forces  in  K-set  (Integer) 

REUNMK  Relation  generated  by  the  UZAERO  module  (Input) 

[PHIKH]  Eigenvectors  in  the  K-set  (Input) 

[QHG(BC,SUB)]  Generalized  gust  forces  (Output) 

Application  Calling  Sequence: 

None 


Method: 

The  module  reads  flutter  identification  number  for  the  subcase  SUB  from  relation  CASE  and 
then  extracts  gust  ID  number  from  relation  ASESOL.  If  it  is  zero  the  module  stops;  otherwise 
matrix  [QHG(BC,SUB)]  is  calculated. 
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Engineering  Application  Module:  MGFL 
Entry  Point:  MGFL 


Purpose: 

To  update  the  generalized  dynamic  and  aerodynamic  matrices,  and  sensitivities  of  the 
generalized  matrices,  for  the  modified  structure. 

MAPOL  Calling  Sequence: 

CALL  MGFL  ( NITER,  BC,  SUB,  HSIZE(BC),  ESIZE(BC),  LAMBDA,  [Mil],  [PSI], 

[GM2(BC,SUB)],  [GK2(BC,SUB)],  [GB2(BC,SUB)],  [MHHFL(BC,SUB)], 
[KRHHFL(BC,SUB)],  [BHHFL(BC,SUB)],  [DTMSB(BC,SUB)], 

[EMSB (B C, SUB)] ,  [AMSB(BC,SUB)],  [DTMS(BC,SUB)],  [EMS(BC,SUB)], 
[AMS (BC, SUB)],  [MIC(BC,SUB)],  [MHC(BC,SUB)],  [DGMVB(BC)], 
[DGKVB(BC)],  [DGMCB(BC)],  [DGMV(BC)],  [DGKV(BC)], 
[DGMVC(BC)],); 


NITER 

Design  iteration  number  (Integer,  Input) 

BC 

Boundary  conditions  identification  number  (Integer,  Input) 

SUB 

Flutter  subcase  number  (Integer,  Input) 

HSIZE(BC) 

The  number  of  Eigenvectors  computed  (Integer,  Input) 

ESIZE(BC) 

The  number  of  extra  point  degrees  of  freedom 
(Integer,  Input) 

LAMBDA 

Relation  of  normal  modes  eigenvalues  output  from  the 
REIG  module  (Input) 

[Mil] 

Modal  mass  matrix  created  in  the  REIG  module  (Input) 

[PSI] 

Matrix  of  Eigenvectors  of  the  reduced  modal  analysis 
(Input) 

[GM2(BC,SUB)] 

Generalized  direct-input  mass  matrix  (Input) 

[GK2(BC,SUB)] 

Generalized  direct-input  stiffness  matrix  (Input) 

[GB2(BC,SUB)] 

Generalized  direct-input  damping  matrix  (Input) 

[MHHFL(B  C,  SUB)] 

Modal  mass  matrix  (Output) 
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[KRHHFL(BC,SUB)] 

Real  modal  flutter  stiffness  matrix  (Output) 

[BHHFL(B  C,  SUB)] 

Modal  flutter  damping  matrix  (Output) 

[DTMSB  (B  C,  SUB)] 

Matrix  [Z)T]  of  aerodynamic  approximation  for  the  baseline 
(Input) 

[EMSB(BC,SUB)] 

Matrix  [E]  of  aerodynamic  approximation  for  the  baseline 
(Input) 

[AMSB(BC,SUB)] 

Matrix  which  consists  of  matrices  [A{\  of  aerodynamic 
approximation  for  the  baseline  (Input) 

[DTMS(BC,SUB)] 

Matrix  [Z)T]  of  aerodynamic  approximation  for  the  modified 
structure  (Output) 

[EMS(BC,SUB)j 

Matrix  [£]  of  aerodynamic  approximation  for  the  modified 
structure  (Output) 

[AMS(BC,SUB)] 

Matrix  which  consists  of  matrices  [A{\  of  aerodynamic 
approximation  for  the  modified  structure  (Output) 

[MIC(BC,SUB)] 

Mass  coupling  matrix  for  structural  and  control  modes  (Input) 

[MHC(B  C,  SUB)] 

Mass  coupling  matrix  for  structural  and  control  modes  in  new 
modal  basis  (Output) 

[DGKVB] 

A  matrix  of  design  sensitivity  of  the  baseline  generalized 
stiffness  matrix  (Input) 

[DGMVB] 

A  matrix  of  design  sensitivity  of  the  baseline  generalized  mass 
matrix  (Input) 

[DGMCB] 

A  matrix  of  design  sensitivity  of  mass  coupling  between 
structural  and  control  modes  of  the  baseline  (Input) 

[DGKV] 

A  matrix  of  design  sensitivity  of  the  modified  generalized 
stiffness  matrix  (Output) 

[DGMV] 

A  matrix  of  design  sensitivity  of  the  modified  generalized  mass 
matrix  (Output) 

[DGMVC] 

A  matrix  of  design  sensitivity  of  mass  coupling  between 
structural  and  control  modes  (Output) 
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Application  Calling  Sequence: 


None 

Method: 

The  module  reads  the  generalized  stiffnesses  from  relation  LAMBDA  and  builds  a  diagonal 
matrix  which  is  the  first  term  in  Eq.  (2.24)  of  the  Theoretical  Manual.  Then  the  module  reads 
the  set  identification  number  of  the  user-input  damping  parameters  (namely  entities  VSDAMP 
and  TABDMP1)  for  the  current  flutter  subcase  from  relation  CASE.  If  it  is  zero,  the  first  and 
third  terms  in  Eq.  (2.26)  of  the  Theoretical  Manual  are  eliminated.  Otherwise  the  user-input 
damping  matrix,  which  is  the  sum  of  these  terms,  is  created.  The  module  reads  the  names  of 
K2PP,  M2PP,  and  B2PP  for  the  current  flutter  subcase  from  relation  CASE.  If  a  name  is  not 
empty,  the  corresponding  term  dealing  with  a  direct-input  matrix  (stiffness,  mass,  or  damping) 
is  calculated.  Finally,  the  updated  matrices  are  created  according  to  Eq.  (2.24)  -  (2.26)  of  the 
Theoretical  Manual.  Then  the  aerodynamic  approximation  matrices  are  updated,  matrices  [Z)T] 
and  [A]  are  post-multiplied  by  [XF],  matrices  [A{\  pre-multiplied  by  [T11]  and  post-multiplied  by 
[T*]  according  to  the  Eq.  (2.42)  of  the  Theoretical  Manual.  Sensitivities  of  the  generalized 
matrices  are  also  pre-multiplied  by  [VPT]  and  post-multiplied  by  ['F]  according  to  Eq.  (2.30), 
(2.3 1)  of  the  Theoretical  Manual. 

Design  Requirements: 

1.  None 

Error  Conditions: 


1.  None. 
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Engineering  Application  Module:  MIST 
Entry  Point:  MIST 


Purpose: 

To  approximate  the  unsteady  aerodynamic  force  coefficient  matrix  [Q(ik)J  by  the  rational 
function  [  Q(ik )]  =  [AJ  +  ik[AJ  -  k2  [  A2J  +  ik[D ](ik[  I]- [R])~'[E]  aerodynamic 
force. 


MAPOL  Calling  Sequence: 

CALL  MIST  (NITER,  BC,  SUB,  GSIZEB,  [QHHLFL(B C, SUB)],  [QHC(BC,SUB)], 
[QHG(BC,SUB)],  [Mil],  [MCC(BC,SUB)],  [MHC(BC,SUB)], 
LAMBDA,  [PHIG(BC)],  [DTMS(BC,SUB)],  [EMS(BC,SUB)], 
[AMS(BC,SUB)],  [RMS(BC,SUB)],  [NLCR(BC,SUB)] ); 


NITER 

BC 

SUB 

GSIZEB 

[QHHLFL(BC,  SUB)] 
[QHC(BC,SUB)] 
[QHG(B  C,  SUB)] 
[Mil] 

[MCC(BC,SUB)] 

[MHC(BC,SUB)] 

LAMBDA 

[PfflG(BC)] 

[DTMS(BC,SUB)] 


Design  iteration  number  (Integer,  Input) 

Boundary  conditions  identification  number  (Integer,  Input) 

Flutter  subcase  number  (Integer,  Input)} 

The  size  of  the  structural  set  (Integer,  Input) 

Generalized  aerodynamic  matrix  (Input) 

Generalized  aerodynamic  control  force  (Input) 

Generalized  aerodynamic  gust  force  (Input) 

Modal  mass  matrix  created  in  the  REIG  module  (Input) 

Mass  coupling  matrix  for  control  modes  (Input) 

Mass  coupling  matrix  between  structural  and  control  modes 
(Input) 

Relation  of  normal  modes  eigenvalues  output  from  the  [REIG] 
module  (Input) 

Eigenvalues  in  the  global  set  (Input) 

Matrix  [ D T]  of  aerodynamic  approximation  (Output) 
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[EMS(BC,SUB)] 

[AMS(BC,SUB)] 

[RMS(BC,SUB)J 

[NLCR(BC,SUB)] 


Matrix  [£]  of  aerodynamic  approximation  (Output) 

Matrix  which  consists  of  matrices  [,4;]  of  aerodynamic 
approximation  (Output) 

Vector  of  lag  terms  of  aerodynamic  approximation  (Output) 
Vector  of  column  ranges  for  Roger  approximation  (Input) 


Application  Calling  Sequence: 
None 


Method: 

The  method  first  extracts  FLCOND  ID  of  the  SUB  flutter  case  from  the  CASE  relation  and 
then  extracts  RAAID,  CONID  and  GUSID  from  the  associated  ASESOL  relation.  General 
approximation  parameters  are  then  extracted  from  the  MINST  AT,  AEROLAG,  DINIT  and 
APCONST  relations.  The  generalized  matrices  are  normalized  to  unit  generalized  masses  of 
the  structural  modes.  Physical  weighting  is  performed  with  the  parameters  of  the  PWEIGHT 
relation  if  APNID  >  0  in  the  MINSTAT  relation.  Approximation  constraints  are  applied 
according  to  the  APCONST  relation.  If  the  ITMAX  parameter  in  MINSTAT  is  zero,  the 
approximation  is  done  by  Roger's  method.  Otherwise,  ITMAX  [Z>]— iterations  are 
performed  to  construct  the  minimum-state  approximation.  The  printout  includes  the  matrix 
[<9*]  of  the  maximal  weighted  aerodynamic  data  terms,  and  the  history  of  weighted 
approximation  error. 

Design  Requirements: 

1.  Can  be  called  only  after  FLUTQHHL. 

Error  Conditions: 


1.  None. 
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Engineering  Application  Module:  RMODE 

Entry  Point:  RMODE 

Purpose: 

To  perform  approximated  normal  modes  analysis. 

MAPOL  Calling  Sequence: 

CALL  RMODE  (  NITER,  BC,  ND,  USET(BC),  [GK(BC)], 

[GM(BC)],  [MRR],  [D],  LAMBDA,  [PSI],  [Mil],  INORM );  V} 


NITER 

BC 

ND 

USET(BC) 

[GK(BC)] 

[GM(BC)] 

[MRR] 

[D] 

LAMBDA 

[PSI] 

[Mil] 

INORM 


Design  iteration  number  (Integer,  Input) 

Boundary  conditions  identification  number  (Integer,  Input) 

Desired  number  of  roots  in  approximated  modal  analysis  (Integer, 
Input) 

The  entity  defining  structural  sets  for  the  current  boundary  conditions 
(Input) 

Generalized  stiffness  matrix  (Input)} 

Generalized  mass  matrix  (Input)} 

The  rigid-body  mass  matrix  (Input)} 

The  rigid-body  transformation  matrix  (Input)} 

Relation  of  normal  modes  eigenvalues  output  from  the  REIG  module 
(Output) 

Matrix  of  Eigenvectors  of  approximated  modal  analysis  (Output) 

Modal  mass  matrix  created  in  the  REIG  module  (Output)} 

Flag  of  mode  normalization  required  in  EIGR  bulk  card:  INORM=l  - 
POINT;  INORM=2  -  MASS;  INORM=3  -  MAX  (Integer,  Output) 


Application  Calling  Sequence: 
None 


25 


Method: 


The  module  reads  from  relation  CASE  the  set  identification  number  of  EIGR  entity  for  the 
current  boundary  condition.  Then  it  opens  this  EIGR  entity,  defines  the  flag  of  mode 
normalization  INORM  and  updates  the  entity  as  follows:  METHOD='MGIV';  MINFREQ=0.0; 
MAXFREQ=0.0;  ROOTEST1=0;  ROOTDESl=ND;  NORM=’MASS';  SETID  and 
ORTHPARM  remain  with  no  change.  Finally  the  module  calls  application  module  REIG  for 
performing  normal  modes  analysis. 

Design  Requirements: 

1.  None. 


Error  Conditions: 


Engineering  Application  Module:  UPDV 

Entry  Point:  UPDV 

Purpose: 

To  update  design  variables. 


MAPOL  Calling  Sequence: 

CALL  UPDV  (  VTNIG,  GLBDES,  NDV,  DVFLG,  [DELDV]  ); 


VINIG 

Relation  of  updated  values  of  global  design  variables  (Input) 

GLBDES 

Relation  of  global  design  variables  (Input) 

NDV 

Number  of  design  variables  (Integer,  Input) 

DVFLG 

Logical  flag  indicating  that  design  variables  were  updated  via  VINIG 
bulk  data  entries  (Logical,  Output) 

[DELDV] 

A  vector  of  differences  between  updated  design  variables  and  those 
used  for  building  modal  data  base  (Output) 

Application  Calling  Sequence: 

None 

Method: 

The  module  retrieves  the  updated  values  of  the  design  variables  from  relation  VINIG.  If  the 
number  of  entries  is  zero,  logical  flag  UPDV  is  set  to  FALSE  and  the  module  quits. 
Otherwise,  logical  flag  UPDV  is  set  to  TRUE,  design  variables  values  are  retrieved  from 
relation  GLBDES,  and  the  vector  of  differences  between  updated  values  and  those  retrieved 

from  relation  GLBDES  is  calculated.  Finally  the  values  of  design  variables  in  relation  GLBDES 
are  replaced  by  the  updated  values. 

Design  Requirements: 

1.  None. 

Error  Conditions: 

1.  None. 
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Chapter  3  Database  Entity  Descriptions 


This  chapter  includes  a  complete  description  of  added  ASE  Data  Base  entities. 

Entity:  AMS 

Entity  Type:  MATRIX 

Description:  Matrix  which  consists  of  matrices  [AO],  [Al],  [A2]  of  aerodynamic 
approximation  for  modified  structure. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  total  number  of  computed  structural, 
control,  and  gust  modes  multiplied  by  3.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 

Created  By:  MIST 

Entity:  AMSB 

Entity  Type:  MATRIX 

Description:  Matrix  which  consists  of  matrices  [AO],  [Al],  [A2]  of  aerodynamic 
approximation  for  baseline. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  total  number  of  computed  structural, 
control,  and  gust  modes  multiplied  by  3.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 

Created  By:  MIST 

Entity:  DGKV 

Entity  Type:  MATRIX 

Description:  Contains  the  design  sensitivities  of  the  generalized  stiffness  matrix  for 
modified  structure. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  number  of  computed  structural  modes 
times  the  number  of  design  variables.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 
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Created  By:  MGFL 


Entity:  DGKVB 

Entity  Type:  MATRIX 

Description:  Contains  the  design  sensitivities  of  the  generalized  stiffness  matrix  for 
baseline. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  number  of  computed  structural  modes 
times  the  number  of  design  variables.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 

Created  By:  MAPOL 

Entity:  DGMV 

Entity  Type:  MATRIX 

Description:  Contains  the  design  sensitivities  of  the  generalized  mass  matrix  for  modified 
structure. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  number  of  computed  structural  modes 
times  the  number  of  design  variables.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 

Created  By:  MGFL 

Entity:  DGMVB 

Entity  Type:  MATRIX 

Description:  Contains  the  design  sensitivities  of  the  generalized  mass  matrix  for  baseline. 

Matrix  Form:  The  number  of  columns  is  equal  to  the  number  of  computed  structural  modes 
times  the  number  of  design  variables.  The  number  of  rows  is  equal  to  the 
number  of  computed  structural  modes. 

Created  By:  MAPOL 
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Entity: 

Entity  Type: 
Description: 

Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 

Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 


DGMCB 

MATRIX 

Contains  the  design  sensitivities  of  the  mass  coupling  matrix  between 
structural  and  control  modes  for  baseline. 

The  number  of  columns  is  equal  to  the  number  ofcontrol  modes  times  the 
number  of  design  variables.  The  number  f  rows  is  equal  to  the  number  of 
computed  structural  modes. 

MAPOL 

DGMVC 

MATRIX 

Contains  the  design  sensitivities  of  the  mass  coupling  matrix  between 
structural  and  control  modes  for  modified  structure. 

The  number  of  columns  is  equal  to  the  number  of  control  modes  times  the 
number  of  design  variables.  The  number  of  rows  is  equal  to  the  number  of 
computed  structural  modes. 

MGFL 

DTMS 

MATRIX 

Matrix  [DT]  of  aerodynamic  approximation  for  modified  structure. 

The  number  of  columns  is  equal  to  the  number  computed  structural  modes. 
The  number  of  rows  is  equal  to  the  number  of  aerodynamic  states. 

MIST 
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Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 


DTMSB 

MATRIX 

Matrix  [DT]  of  aerodynamic  approximation  for  baseline. 

The  number  of  columns  is  equal  to  the  number  of  computed  structural  modes. 
The  number  of  rows  is  equal  to  the  number  of  aerodynamic  states. 

:MIST 


EMS 

MATRIX 

Matrix  [E]  of  aerodynamic  approximation  for  modified  structure. 

The  number  of  columns  is  equal  to  the  total  number  of  computed  structural, 
control,  and  gust  modes.  The  number  of  rows  is  equal  to  he  number  of 
aerodynamic  states. 

MIST 

EMSB 

MATRIX 

Matrix  [E]  of  aerodynamic  approximation  for  baseline. 

The  number  of  columns  is  equal  to  the  total  number  of  computed  structural, 
control,  and  gust  modes.  The  number  of  rows  is  equal  to  the  number  of 
aerodynamic  states. 

MIST 

GB2 

MATRIX 

Contains  the  generalized  direct-input  damping  matrix. 
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Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 


Real  matrix  having  one  row  and  one  column  for  each  computed  structural 
mode. 

Module  BDMAT 

GK 

MATRIX 

Contains  the  generalized  stiffness  matrix  for  modified  structure. 

Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
Module  GMGKUPD 

GK2 

MATRIX 

Contains  the  generalized  direct-input  stiffness  matrix. 

Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
Module  BDMAT 

GKB 

MATRIX 

Contains  the  generalized  stiffness  matrix  for  baseline  structure. 

Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
Module  GMAT 

GM 

MATRIX 

Contains  the  generalized  mass  matrix  for  modified  structure. 
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Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 
Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 


Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
Module  GMGKUPD 

GM2 

MATRIX 

Contains  the  generalized  direct-input  mass  matrix. 

Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
Module  BDMAT 

GMB 

MATRIX 

Contains  the  generalized  mass  matrix  for  baseline  structure. 

Real  matrix  having  one  row  and  one  column  or  each  computed  structural  mode. 
Module  GMAT 

KRHHFL 

MATRIX 

Real  stiffness  matrix  in  the  modal  set  used  in  ASE  flutter  analysis. 

Real  matrix  having  one  row  and  one  column  for  each  computed  structural  mode. 
BDMAT,  MGFL 

MCC 

MATRIX 

Mass  matrix  for  control  modes. 

Real  matrix  having  one  row  and  one  column  for  each  control  mode. 
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Created  By: 

Entity: 

Entity  Type: 
Description: 

Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 

Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 

Matrix  Form: 

Created  By: 


CNTMOD 

MIC 

MATRIX 

Matrix  of  mass  coupling  between  structural  and  control  modes  of  the  modified 
structure  in  the  modal  basis  of  the  baseline. 

Real  matrix  having  one  column  for  each  control  mode  and  one  row  for  each 
computed  structural  mode. 

GMGKUPD 

MICE 

MATRIX 

Matrix  of  mass  coupling  between  structural  and  control  modes  of  the  baseline 
structure. 

Real  matrix  having  one  column  for  each  control  mode  and  one  row  for  each 
computed  structural  mode. 

CNTMOD 

MHC 

MATRIX 

Matrix  of  mass  coupling  between  structural  and  control  modes  of  the  modified 
structure  in  the  modified  modal  basis. 

Real  matrix  having  one  column  for  each  control  mode  and  one  row  for  each 
computed  structural  mode. 

MGFL 
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Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 


NLCR 

MATRIX 

Contains  the  column  ranges  for  Roger  approximation. 

The  number  of  columns  is  equal  to  1.  The  number  of  rows  is  equal  to  number 
of  computed  structural  modes  plus  one. 

MIST 

PHFLO 

MATRIX 

Contains  the  numbers  of  structural  modes  to  be  used  in  ASE  analysis. 

The  number  of  columns  is  equal  to  1.  The  number  of  rows  is  equal  to  number 
of  structural  modes  to  be  used  in  ASE  analysis. 

ASEFLUT 

PHIAB 

MATRIX 

Contains  the  Eigenvectors  of  the  baseline  modes  in  the  a-set  set. 

Real  rectangular  matrix  having  one  column  for  each  computed  structural  mode 
and  one  row  for  each  a-set  degree  of  freedom  . 

MAPOL 

PLROG 

MATRIX 

Truncation  vector  for  Roger  approximation. 

The  number  of  columns  is  equal  to  1.  The  number  of  rows  is  equal  to  number 
of  structural  modes  to  be  used  in  ASE  analysis  times  number  of  aerodynamic 
states. 
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Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 

Entity: 

Entity  Type: 
Description: 
Matrix  Form: 

Created  By: 


ASEFLUT 

PSI 

MATRIX 

Contains  the  Eigenvectors  of  the  reduced  modal  analysis. 

Real  rectangular  matrix  having  one  column  and  one  row  for  each  computed 
structural  mode. 

RMODE 

QHC 

MATRIX 

Unsteady  generalized  aerodynamic  control  forces. 

Complex  matrix  with  number  of  columns  being  equal  to  the  product  of  the 
number  of  control  surfaces  and  the  number  of  M-k  pairs  for  which 
aerodynamics  are  required  and  number  of  rows  being  equal  to  the  number  of 
computed  structural  modes. 

CNTMOD 

QHG 

MATRIX 

Unsteady  generalized  aerodynamic  gust  forces. 

Complex  matrix  with  number  of  columns  being  equal  to  the  productof  the 
number  of  gust  modes  and  the  number  of  M-k  pairs  for  which  aerodynamics 
are  required  and  number  of  rows  being  equal  to  the  number  of  computed 
structural  modes. 

GSTMOD 
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Entity:  RMS 

Entity  Type:  MATRIX 

Description:  Vector  of  the  aerodynamic  lag  terms. 

Matrix  Form:  The  number  of  columns  is  equal  to  1.  The  number  of  rows  is  equal  to  the 
number  of  aerodynamic  states. 

Created  By:  MIST 

Entity:  ASELMBD 

Entity  Type:  RELATION 

Description:  Contains  results  of  an  ASE  flutter  analysis  for  a  series  of  boundary 
conditions,  Mach  numbers  and  atmospheric  densities. 

Relation  Attributes: 


NAME 

TYPE/KEY 

DESCRIPTION 

NITER 

Integer 

Iteration  number 

BCID 

Integer 

The  boundary  condition  number 

MACH 

Real 

Mach  number  of  the  flutter  analysis 

RHOREF 

Real 

Reference  atmospheric  density 

RHO 

Real 

Atmospheric  density  of  the  flutter  analysis 

VELOCITY 

Real 

True  velocity  of  the  flutter  analysis 

FSID 

Integer 

Flutter  set  identification 

SCNUM 

Integer 

Flutter  subcase  identification  number 

MODENO 

Integer 

Mode  number  associated  with  flutter 

RLAMB 

Real 

Real  part  of  the  flutter  eigenvalue 

ILAMB 

Real 

Imaginary  part  of  the  flutter  eigenvalue 

DAMPVAL 

Real 

Damping  ratio 

OMEGA 

Real 

Frequency  in  radians  per  second  of  the  flutter 
eigenvalue 

=2*VELOCITY*ILAMB/REFB 

PNUM 

Integer 

Pointer  to  CONST  tuple  for  the  associated 
constraint 

Created  By:  Module  ASEFLUT 
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Entity: 


ASEMARG 


Entity  Type:  RELATION 

Description:  Contains  results  of  an  ASE  SISO  margin  analysis,  ASE  MIMO  margin 
analysis,  and  open-loop  vehicle  gain  calculation  for  a  series  of  boundary 
conditions,  Mach  numbers  and  atmospheric  densities. 

Relation  Attributes: 


NAME 

TYPE/KEY 

DESCRIPTION 

NITER 

Integer 

Iteration  number 

BCID 

Integer 

The  boundary  condition  number 

SCNUM 

Integer 

Flutter  subcase  identification  number 

FSID 

Integer 

Flutter  set  identification 

MACH 

Real 

Mach  number  of  the  flutter  analysis 

RHOREF 

Real 

Reference  atmospheric  density 

DENRAT 

Real 

Density  ratio  of  the  flutter  analysis 

VELOCITY 

Real 

True  velocity  of  the  flutter  analysis 

GROW 

Integer 

Row  number  for  the  gain  of  the  SISO  margin 
analysis 

GCOL 

Integer 

Column  number  for  the  gain  of  the  SISO 
margin  analysis 

STAB 

Integer 

Flag  indicating  stability  of  the  closed-loop 
system.  See  Remark  1 

MARTYPE 

Text  (4) 

Margin/constraint  type.  See  Remark  2 

MARINF 

Integer 

Flag  indicating  infinity  value  of  SISO  margin. 
See  Remark  3 

VALDIM 

Real 

Margin  or  open-loop  vehicle  gain  in  relative 
units 

MARVAL 

Real 

Value  in  dB  or  degrees  (for  SISO  margins). 
See  Remark  4 

OMEGAC 

Real 

Cross-over  frequency  in  rad/s.  See  Remark  5 

CONVAL 

Real 

Constraint  value 

PNUM 

Integer 

Pointer  to  CONST  tuple  for  the  associated 
constraint 

Created  By:  Module  ASEFLUT 
Remarks: 

1 .  The  stability  flag  is: 

1  =  Closed-loop  system  is  stable 
-1  =  Closed-loop  system  is  unstable 
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2.  The  margin/constraint  types  are: 

UGM  =  Upper  gain  margin 
LGM  =  Lower  gain  margin 
UPM  =  Upper  phase  margin 
LPM  =  Lower  phase  margin 

S  VI  =  Singular  value  for  input 
SVO  =  Singular  value  for  output 
GC  =  Gain  constraint 

3.  The  infinity  flag  is: 

0  =  SISO  margin  is  finite 
1  =  SISO  margin  is  infinite 

4.  For  margin  types  SVI,  SVO,  GC  attribute  MARVAL  contains  the  same  value  as  VALDIM. 

5.  For  MIMO  margins  attribute  OMEGAC  contains  frequency  for  which  the  singular  value  was 
calculated. 
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